SYMPLECTIC STRUCTURES FOR THE CUBIC SCHRODINGER 
EQUATION IN THE PERIODIC AND SCATTERING CASE. 



K.L. VANINSKY 

Abstract. We develop a unified approach for construction of symplectic forms 
for ID integrable equations with the periodic and rapidly decaying initial data. 
As an example we consider the cubic nonlinear Schrodinger equation. 



1. Introduction 

1.1. General remarks. The main technical tool for the study of soliton systems is 
commutator formalism. All fashionable soliton systems like the Korteveg-de-Vriez 
equation (KdV), the cubic nonlinear Schrodinger equation (NLS), the sin-Gordon 
equation, the Toda lattice, etc; have such representation. Within the commutator 
formalism approach the dynamical system appears as a compatibility condition 
for an over-determined system of equations. As an example, we consider the NLS 
equation with repulsive nonlinearity 1 

= _^" + 2|?/>| 2 7/s 

where tp(x,t) is a complex function of spatial variable x and time t. The flow is a 
compatibility condition for the commutator 



with 



and 



[d t -V 3 ,d x -V 2 ] = 0, 



xr / ,\ T/ *A iX ( 1 \ _ / V 

V 2 (x,t)=V = --a 3 + Y = - Y l _ x j + f I 



A 2 

V 3 (x,t) = - ia 3 - XY + \4>\ 2 i<J 3 - 10^. 



The corresponding auxiliary linear problem 

(d x - V)f = 0, / 



h 
h 
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can be written in the form of an eigenvalue problem for the Dirac operator 



Another important feature of soliton systems is the Hamiltonian formulation. 
Here we assume that the potential tp(x,t) is 2/-periodic: ip(x + 21, t) = ip(x,t). 
For instance, the NLS flow can be written as 



The bracket is nondegenerate. The corresponding symplectic form (up to a scalar) 
is: 



A priori it is not clear why the dynamical system, which arises as a compatibility 
condition has a Hamiltonian formulation. To put it differently, is it possible to 
obtain Hamiltonian formalism from the spectral problem? 

Here we would like to make some historical remarks. Originally, the Hamiltonian 
formulation of basic integrable models was found as an experimental fact. For the 
KdV equation the computation of symplectic structure in terms of the scattering 
data was performed by Faddeev and Zakharov, [5]. It involved some nontrivial 
identities for the products of solutions. Later Kulish and Reiman, [14], noted 
that all higher symplectic structures also can be written in terms of the scattering 
data. Again, they used the scheme of [5] and explicit calculations. Finally, we 
note that Zakharov and Manakov, [28], for the NLS equation adopted a different 
approach. Instead of the symplectic structure they worked with the corresponding 
Poisson bracket. Again, using explicit formulas for the product of solutions they 
computed the Poisson bracket between the coefficients of the scattering matrix. An 
appearance of explicit formulas that are the moving force of all these computations 
seems to be quite mysterious. This was already discussed in the literature [4], and 
described as a "computational miracle". 

The standard assumption needed to carry out spectral analysis is that the po- 
tential either is periodic or has rapid decay at infinity. We refer to the latter case 
as scattering. Recently, in connection with the Seiberg-Witten theory [23, 24], 
Krichever and Phong, [13] developed a new approach for the construction of sym- 
plectic formalism. The latest exposition of their results can be found in [3]. The 
main idea of the Krichever-Phong approach is to introduce in a universal way the 
two-form on the space of auxiliary linear operators. This form is written in terms 




with Hamiltonian H 3 — | J l _ t \ip'\ 2 + \ip\ 4 dx =energy and bracket 
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of the operator itself and its eigenfunctions. The goal of this paper is to review 
the Krichever-Phong approach in the case of ID periodic NLS and to extend it to 
the scattering case. Within the unified approach, we reduce the number of formu- 
las and eliminate unnecessary explicit computations. For instance computation of 
the symplectic form in terms of the spectral data (both in the periodic and the 
scattering case) becomes an application of the Cauchy residue theorem. 

1.2. The periodic case. We assume that the potential is periodic with the period 
21 : ip(x + 21, t) = ip(x,t). The Krichever-Phong formula, in the NLS context, 
takes the form 



This formula defines a closed 2-form u on the space of operators d x — V(x, A) 
with 21 periodic potential. The set-up for this formula is broadly as follows 2 . 
The eigenvalue problem 



has special solutions, so-called Floquet solutions determined by the property e(x+ 
21, A) = w(X)e(x, A). The complex constant w(X) is called a Floquet multiplier. For 
each value of the spectral parameter A there are two linearly independent Floquet 
solutions and two distinct Floquet multipliers. These solutions and correspond- 
ingly multipliers become single-valued functions of a point on the two-sheeted 
covering of the plane of spectral parameter A. The simple points of the peri- 
odic/antiperiodic spectrum of the eigenvalue problem constitute branching points 
of the cover. We assume that there is a finite number of simple points (so-called, 
finite gap potential). 

This two sheeted covering constitutes a hyperelliptic Riemann surface T with two 
infinities P + and P_ (Fig. 1). Each point Q = (A, ±) of T is specified by the value 
of spectral parameter A and the sheet "+" or "-" which determines the Floquet 
multiplier w(Q) corresponding to this A. At every point of the curve we also have 
a Floquet solution e(x, Q) which becomes a function of the point Q and satisfies 
the identity e[x + 21, Q) = w(Q)e(x,Q). The Floquet solution e(x, Q) has an 
exponential singularity at infinities and plays the role of so-called Baker-Akhiezer 
function for the curve V . 

At every point of the curve V we can define another solution e*(x, Q). This is the 
Floquet solution which is brought from a point on the different sheet but with the 
same value of the spectral parameter A. It is transposed and suitably normalized. 
The operator Jd x — JV(x,\) acts on the solution e*(x, Q) as an adjoint i.e. on 




res < e*J 5V A 5e > dX. 



i'± 





2 We refer to Section 2 for detailed discussion. 
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Figure 1. Smooth Riemann surface T. 



the right: 

e*(x,Q) \Jd x - JV] = 0. 

It is assumed that the phase space consists of smooth 2/-periodic functions 
ip(x) or equivalently operators d x — V(x, A) with 2/-periodic potential. The NLS 
flow acts on this space as it acts on the space of functions ip(x). All notions 
of differential geometry with obvious conventions can be applied to this space of 
operators. On the space of potentials we have a variation 5ip(x). Thus for a 
fixed value of the spectral parameter A we have well defined variation 5V(x, A). 
The variations Se(x, Q), Se*(x, Q) are defined correctly when A = A(Q) is fixed. 
Therefore, at each point Q of the surface T we have well defined meromorphic in 
Q the two-form 

< e*J5V A 5e > d\. 

It takes values in the space of skew-symmetric two-forms on the space of operators 
d x — V. The result of Krichever and Phong states that the sum of residues of this 
form at infinities P± is nothing but the symplectic form ujq. 

The formula has a lot of good properties. First, it produces all higher symplectic 
structures by introducing the weight A n under the residue 

u n = ^2 res X 1 < e *J SV A 5e > d X, n = 1, 2, 

p± 

Second, it easily leads to the Darboux coordinates, or in physics terminology the 
separation variables, see Sklyanin [25]. These are local coordinates where the 
symplectic form ojq takes the simple canonical form 

2 . 

= - 8p(lk) A <5A(7 fe ). 
k 
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This merits special explanation. It is well-known since the work of Flashka- 
McLaughlin, [6], that the poles 7 fc of Floquet solutions lead to the Darboux co- 
ordinates for symplectic forms 3 . Recently, a lot of work was performed, [15], to 
construct such variables for the Ruijsenaars-Sneider and the Moser-Calogero sys- 
tems. This required formidable technical machinery and extensive computations. 
At the same time, as it was demonstrated by Krichever, [11], the formula leads to 
the same result only by applying the Cauchy residue theorem. 

1.3. The Scattering case. The main goal of the present paper is to show that 
suitably interpreted the new approach can be adopted for soliton systems with 
rapidly decaying initial data on the entire line. This is so-called the scattering 
case 4 . 

For such potentials one can define so-called Jost solutions J±(x, A). These are 
matrix solutions of the auxiliary linear problem J' ± = V J± with the asymptotics 

J'±(x, A) = exp (— i|xcr 3 ) + o(l), as x — > ±oo. Their columns J± 

analytic in the corresponding upper/lower half-plane. 



-•(i) ~(2) 



are 




Figure 2. Singular Riemann surface T, 



Our construction of the associated Riemann surface is a geometrical interpre- 
tation of what is called the Riemann-Hilbert approach to the scattering problem, 
see [4] . A singular curve is obtained by taking two copies of the complex plane 
and gluing them to each other along the real line (Fig. 2). The curve has two 
infinities P + and P_ and continuum set of singular points above the real line. The 

3 See also Novikov-Veselov, [22], for general discussion. 
4 We refer to Section 3 for detailed definitions. 
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standard Jost solutions are lifted on and become the single valued function of 
a point on the curve. Different branches of BA function are connected along the 
real line by the scattering matrix S: 



S(X) 



1 b 
-b 1 



The Jost solution has exponential singularity at infinities and plays the role of the 
Baker- Akhiezer function for the curve T^. This construction is explained in detail 
in Section 3. 

The formula of Theorem 3.4 looks similar to the periodic case 

u = trace res ]- [< H* + J5V A 5H + > + < H*_J8V A 5H_ >] dX. 
The only difference now is that we work with the matrix solutions 



H + (X)= j { l\X),j?(X) 



;(2)» 



and 



H-(\)= j+ > (X),j [ L , (X) 



;( 2 )/ 



and H+(X) = a t Hl, H±(\) = ^Hl, with 



0"! 



1 

1 



The averaging now corresponds to the integration on the entire line 



< • > 



/+oo 
dx. 
-oo 



The residue can be computed explicitly u = 2i < 5ip A 5ip >. Theorem 3.6 states 
that the symplectic structure can be put in the Darboux form 



-I 

m J 



8b(X) A 8b(X) 
l«(A)| 2 ' 



dX, 



where a and b are coefficients of the scattering matrix S. Again identically to 
the periodic case this result is obtained by applying the Cauchy residue theorem. 
Only now the sum of the residues in the affine part of the curve is replaced by its' 
continuous analog. This is the integral which stays in the right hand side of the 
formula. 

The unified approach to construction of symplectic forms produces an interesting 
problem. As we see, the symplectic form constructed in the periodic case has 
two systems of Darboux coordinates. One system is associated with poles of the 
Floquet solution. It is the divisor-quasimomentum Darboux coordinates. Another 
system of Darboux coordinates is the action-angle variables. At the same time in 
the scattering case we know only one system of Darboux coordinates. These are 
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action-angle variables. What is the correct analog of the divisor-quasimomentum 
in the scattering case? This is a subject of future publication, [27]. 

We conclude this introduction by expressing thanks to A. Its, H. McKean and I. 
Krichever for stimulating discussions. We are also greatful to anonymous referee 
for remarks that helped to improve the presentation. 



2. The Periodic Case. 

2.1. The Direct Spectral Problem. We provide here information needed in 
the next section for construction of symplectic forms. We refer to classical books 
[20, 21] for standard facts of spectral theory and algebraic-geometrical approach 
to solitons. 

The NLS equation 

zifj* = -^j" + 2\ij\ 2 ij, (2.1) 

where tp(x,t) is a smooth complex function 2/-periodic in x, is a Hamiltonian 

system 

with the Hamiltonian H = | J l _ t \ip'\ 2 + \tp\ A dx =energy and the bracket 



OA dB OA dB , 

= — ax. 



dtp(x) d^(x) dtp(x) dtP(x) 

The NLS Hamiltonian H = H 3 is one in the infinite series of conserved integrals 
of motion. 



H x = -J \^\ 2 dx, 
1 f l - 



H 3 = \ijjf + \ijj\ 4 dx, etc. 



These Hamiltonians produce an infinite hierarchy of flows e tXm , m = 1,2, 

The first in the hierarchy is the phase flow e tXl generated by the vector field 

X x : V" = {^,-ffi} = -#- 

The phase flow is a compatibility condition for 

[d t -V 1 ,d x -V 2 ] = 0, (2.2) 
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with 5 Vi = f a 3 and 



2 JU V° f/V^O 

We often omit the subscript V — V 2 . The second, translation flow e iX2 generated 
by 

X 2 : V = {$,H 2 }=il)' 

is equivalent to 2.2 with V\ replaced by V 2 - Finally, the third, original NLS flow 
2.1 is a compatibility condition for 2.2 with V\ replaced by 

A 2 

V 3 = —ioz - XY + |Vf^ 3 - i<7 3 *o- 
All flows of infinite hierarhy e tXm , m — 1, 2, . . . commute with each other 

The first times Ti, r 2 and r 3 correspond to the first three flows. 

We introduce a 2 x 2 transition matrix M(x,y, A), x > y; that satisfies 

M'(x, y, A) = V(x, X)M(x, y, A), M(y, y, A) = /. 

The solution is given by the formula 



M(x,y,\)=exp /"V(£,AK. 

J y 



The matrix M(x, y, A) is unimodular because V is traceless. 
The symmetry 

a 1 V(x : A)(Ji = V(x,A) 
produces the same relation for the transition matrix 

a 1 M(x,y,\)(T 1 = M(x,y,\). (2.3) 

Another symmetry 

\/ T (x,A)J= -JV(x,X), 

where J = ia 2) implies 

M T (x,y,X)- 1 J = JM(x,y,X). (2.4) 

The quantity A(A) = |traceM(i, — I, X) is called a discriminant. The formula 
2.3 implies A(A) = A(A) and A(A) is real for real A. The eigenvalues of the 



5 Here and below a denotes the Pauli matrices 



1 \ / -i \ (10 



ai = 1 1 ) ' CT2 = [ i J ' as = V -1 
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monodromy matrix have a name of Floquet multipliers and they are roots of the 
quadratic equation 

w 2 - 2Aw + 1 = 0. (2.5) 

The Floquet multipliers are given by the formula w = A ± \/ A 2 — 1. The values 
of A : w(X) = ±1 constitute the points of the periodic/antiperiodic spectrum. The 
corresponding auxiliary linear problem 

(d x -v)f = o, / r = (/i,/2); 

can be written in the form of an eigenvalue problem for the self-adjoint Dirac 
operator 



i o V / o 

— 1 / \ iip 



The self-adjointness implies that points of the spectra are real. 
Example. Let ip = 0. The corresponding monodromy matrix can be easily com- 
puted M(x,y,X) = e~ l 2 cz ( x ~ v \ We have A(A) = cos XI and double eigenvalues at 
the points A^ = E P L . If n is even/odd, then the corresponding A^ belongs to the 
periodic/ anti-periodic spectrum. 

For a generic potential the double points A^ of the periodic/anti-periodic spec- 
trum split, but they always stay real. The size of the spectral gap is determined, 
roughly speaking, by the corresponding Fourier coefficients of the potential. In our 
considerations we assume that there is a finite number of g + 1 open gaps in the 
spectrum 

• • • < = Ki-l < < A^ < . . . < A n+3 < \^ +g < A n+3+1 = X-n+g+l < ■ ■ ■ 

These are so-called finite gap potentials which are dense among all potentials. 
The Floquet multipliers become single-valued on the Riemann surface: 

r = {Q = (A, w) G C 2 : R(X, w) = det [M(l, -/, A) - wl] = 0}. 

The Riemann surface consists of two sheets covering the plane of the spectral 
parameter A. 

Example. Let ip = 0. We have A(A) = cos XI and quadratic equation 2.5 has the 
solutions w(X) = e ±llx . The Riemann surface V = T + + r_ is reducible and consists 
of two copies of the complex plane C that intersect each other at the points of the 
double spectrum A^. Each part T + or T_ contains the corresponding infinity P + 
or P_. The Floquet multipliers are single valued on T: 

w(Q) = e +iAZ , QeT + ; 
w(Q) = e~ lXl , Qer_. 

For a finite-gap potential the Riemann surface T is irreducible. There are three 
types of important points on T. These are the singular points, the points above 
A = oo and the branch points which we discuss now in detail. 
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• The singular points are determined by the condition 

d x R{X,w) = d w R{X,w) = 0. 

These are the points (X ± , ±1) of the double spectrum. At these points two 
sheets of the curve intersect. 

• There are two nonsingular points P + and P_ above A = oo. At these points 6 

W (Q) = e + ^(l + 0(l/A)), Qe(P + ); (2.6) 

w(Q) = e- ai (l + 0(l/\)), Qe(P-). (2.7) 

• The branch points are specified by the condition 

d w R{X,w) = 0. 

They are different from the singular points and correspond to the simple pe- 
riodic/antiperiodic spectrum. We denote these points by sf = (Xf, (— l) fe ), 
k — n, . . . ,n + g. There are 2(g + 1) of them, each has a ramification index 
2. 

The desingularized curve V is biholomorphicaly equivalent to a hyperelliptic curve 
with branch points at the points of the simple spectrum. We also denote the 
hyperellitic curve by T. The Riemann-Hurwitz formula for the genus of T implies 

R 

genus = — — n + 1, 

where R is a total ramification index and n is the number of sheets. Each branch 
point has a ramification index 1 and therefore R = 2(g + 1) and n — 2. Therefore, 
the genus of T is g, one off the number of open gaps in the spectrum. 
Let e± be a holomorphic involution on the curve T permuting sheets 

e± : (X,w) — > (A, l/w). 

The fixed points of e± are the branch points of Y. The involution e± permutes 
infinities e± : P_ — > P + . Let us also define on V an antiholomorphic involution 

e a : (X,w) — > (X,w). 

The involution e a also permutes infinities and commutes with e±. Points of the 
curve above gaps [A~, A+] where |A(A)| > 1 form g + 1 fixed "real" ovals of e a . We 
call them a-periods. 

The quasimomentum p(Q) is a multivalued function on the curve T. It is intro- 
duced by the formula w(Q) = e tp ^ 21 . Evidently, it is defined up to 2 y i , where n is 
an integer. The asymptotic expansion for p(Q) at infinities can be easily computed 

±p(A) = ^-p ± -y-g---, Qe(P ± ), X — X(Q); 

6 Thc notation Q e (P) means that the point Q is in the vicinity of the point P. 

10 



where p^ = ^j^, k± is an integer and 

Pi = y#i, P2 = jH 2 , p 3 = jH 3 , etc. 

Moreover, the function w(Q) + w(e±Q) does not depend on the sheet and is equal 
to 2A(A). Thus A(A(Q)) = coship(Q)2l and the formula 

1 , 1 A* (A) dX 

implies that differential dp is of the second kind with double poles at the infinities: 
±dp = rf(A + o(l)). The same formula implies that the differential dp is pure 
complex on the real ovals. At the same time, the condition w(s^) = w(s^) requires 
the increment p(s%) — p(s^) to be real. Therefore, dp has zero a-periods 

dp = 0. 

Since the Floquet multiplies are single-valued on V for the 6-periods we have 

nn bk 



/ dp 



I 



n b , eZ, k = 1, 



,9- 



(2.8) 



These are so-called the periodicity conditions, [21]. 
The Floquet solution is the vector-function 

e\x,Q) 



e(x,Q) 



e\x,Q) 



which is a solution of the auxiliary spectral problem e! = Ve with the property 

e(x + 21, Q) = M{1, -/, X)e(x, Q) = w(Q)e(x, Q) (2.9) 
and normalized by the condition 

e 1 (-l,Q) + e 2 (-l,Q) = 1. (2.10) 
Remark. If f(x, A) is a solution of the auxiliary problem 

(d x -V(x,X))f = 

corresponding to A, then / = a±f is a solution of (d x —V(x, A))/ = corresponding 
to A. 

Example. Let ip — 0. The Floquet solution is given by the formula 



e(x,Q)=e +i ^ x+ »e = e +i ^ x+ » 



e(x,Q) 




1 

1 





Qer + , 
Q e r_. 



It has no poles in the affine part of the curve. 
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For a general finite gap potential the situation is more complicated. 
Lemma 2.1. The Floquet solution satisfies the identity 



e(x,e a Q) = a ie {x,Q). 
The Floquet solution e(x, Q) has poles common for both components at the points 

71,72, •••,7s+i ■ 

Projections of poles ji^ = X{jk) are real. Each 7^ lies on the real oval above the 
corresponding open gap [A^T , AjJ"]. Each component e l (x, Q) has g + 1 zeros 

a[(x), a l 2 (x),...,ai +1 (x); i = 1,2. 

These zeros depend on the parameter x. In the vicinity of infinities the function 
e(x, Q) has the asymptotics 



e(x,Q) = e ±i ^ [co/eo + o(l)] , 



Before proceeding to the proof of the Lemma we note that the differential equa- 
tion for the monodromy matrix 



M'{x,y,\) 



-—a 3 + Y 



M{x,y,X), 



M(y,y,X) 



multiplied (gauged) on the left and right by the matrices 



C 



transforms into 

M R (x,y,X) 

with 



1 1 

i —i 



and 



c- 1 



—1 

i 



M R (x,y,X), 



v K — 

^0 — 



M R (y,y,X) = I; 
ip — q + ip. 



q p 

p -q 

This is a real version of the eigenvalue problem which is more convenient in some 
situations, see [19]. The Floquet solution e R (x, Q) corresponding to the real version 
of the eigenvalue problem is related to e(x, Q) by the formula 

e R (x,Q) =Ce(x,Q). 

Therefore the result of the Lemma for e(x, Q) follows from the corresponding 
result for e R (x, Q) given in [19]. We prefer to give a direct proof, though the gauge 
transformation is behind all arguments. 

Proof. The proof is based on the explicit formula for the Floquet solution. Let 



M(x,-l,X) 



12 



(x,-l,X); 



and Mn, M i2 , efc. be the elements of the matrix M(l, —I, A). The Floquet solution 
e(x, Q) is given by the formula 



e(x,Q) 
where A = 
A(Q) 



A(Q) 



m 21 



(x, -l,\) + (i-A(Q)) 



m 12 
m 22 



[x, 



-I A), (2.11) 



X(Q) and the coefficient A(Q) is 
Mi 2 



or 



A(Q) 



w{Q) - M 22 



(2.12) 



M 12 -M u + w(Q) M 2 i - M 22 + 

To prove the formula note that the Floquet solution is a linear combination of 
columns of the monodromy matrix M(x, —I, A): 



e(x,Q)=A(Q) 



m 21 



\x, 



-l,X)+A'(Q) 



m 12 
m 22 



\x, 



-l,X), A = X(Q). 



The normalization condition 2.10 implies A'(Q) = 1 — A(Q). At the same time 
the Floquet solution is an eigenvector of the monodromy matrix 



M(l,-l,X) 



A(Q) 
1 - A{Q) 



w(Q) 



A(Q) 
1-A(Q) 



This leads to two equations 

M U A(Q) + M 12 (l - A(Q)) = w(Q)A(Q), 

or 

M 21 A(Q) + M 22 (l - A{Q)) = w(Q)(l - A{Q)). 

Each equation implies the corresponding formula for A(Q). 
The formulas 2.3 and 2.12 imply 

1 - A(e a Q) = A(Q). 

This and 2.3, 2.11 imply the stated identity for the Floquet solution. 
The relation M R = CMC' 1 implies 



M 



Mf x M 



12 



M« Ml 



'21 1V± 22 

Mn + Mi 2 + M 2 i + M 22 i(M 12 + M 22 - Mn - M 2 i) 
i{M u + M12 - M 2 i - M 22 ) Mn + M 22 - M12 - M 2 i 



Due to 2.3 M R (A) is real for real A. Consider the function Mf 2 (A) and look at 
the roots fx n : Mf 2 (/i n ) =0. For ^ = we have Mf 2 (A) = - sin f with roots at 
the points /i n = n G Z. When we add the potential the roots fi n move but 
stay real. They are caught by open gaps or match double periodic/antiperiodic 
spectrum. Indeed at fi n the matrix M R is lower triangular and real entries M* and 
Ml 



[22 coincide with Floquet multipliers. Since Mj^M^/^) = 1 we have \A(fx n )\ 
Mft + M&l > 1. 
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The points of the divisor ^ k e T lie above the points /i^ on the sheet with 
w(Q) = Mf 2 (fi n ). At these points the denominator in 2.12 

M12 - Mn + w(Q) or M 2 i - M 22 + w(Q) 

vanishes. Indeed from M^(// n ) = we have 

Mia - M n (/i n ) = M 2 i - M 22 (/i n ) 

and iu(Q) = M 2 R 2 (/i„) = M n - M 12 (/i n ) = M 22 - M 21 (/i n ). Moreover M 22 ^ 
M|| These produce a pole of the Floquet solution when /i„ lies in the 
open gap. When /i„ is caught by the periodic/antiperiodic spectrum the ma- 
trix M(l, — i,/i„) = ±1 and the zero of denominator is annihilated by the zero of 
numerator in 2.12. 

The asimptotics of the Floquet solution follows from the formula 2.11 and 

M(x,y,X) = e - l ^ x - y) + o{l), when A -> 00. 

□ 

The Floquet solution e(x, Q) near infinities can be expanded into the asymptotic 



series 



e(x, Q) = e +i ^ x+ » J2 e s(x)\- s = ^ 



s=0 

00 



s=0 

00 



e(x, Q) = e-'f (*+') J2 e s WA- s = {x+l) ^ ^ s A" 
and 6 = 0, ofo = 1- The coefficients b s , d s can be computed from the relation 



d s 
d. 



Qe(P+), 
Qe(P-), 



K 

d' 



b s 



due to the diagonal form of the matrix 



i 




0,1,..., 



(2-13) 



Indeed, relation 2.13 leads to the identities 

-b' s + i)d s = ib s+1 , 

and 

— d' s + ipb s = 0. 
These are supplemented by the boundary condition 

b s (x) + d s (x)\ x= -i = 0, s > 1. 

When s = using 6 = 0, rf = 1 from 2.14 we obtain 

b\ = —iip(x). 
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(2-14) 
(2.15) 
(2.16) 



The identities 2.15 and 2.16 imply 



Similar, we compute 



ip + tp tp 



X 



d 2 



-t' -t 2 



-i L 



X 

o-l^l 2 / 



These formulae will be used for explicit computation of symplectic forms. 
The Floquet solution e(x, Q) satisfies the identity 

[Jd x — JV] e(x, Q) = 0, J = ia 2 ; 

which is just another way to write the spectral problem. Let us define the dual 
Floquet solution e + (x, Q) = [e 1+ (x, Q), e 2+ (x, Q)\ at the point Q as 

e + (x, Q) = e(x,e±Q) T . 

It can be verified by a direct computation that the dual Floquet solution e + (x, Q) 
satisfies 7 

e + (x,Q) [Jd x - JV] = 0. 

The fact that the Wronskian e + (x, Q)Je(x, Q) does not depend on x can be verified 
by differentiation. Introducing the function 



$(Q) = e + (x, Q)Je(x, Q) =< e + (x, Q)Je(x, Q) >= — J e + (x, Q)Je(x, Q)dx, 

we define another dual Floquet solution e*(x,Q) by the formula 

e+(x,Q) 



e*(x,Q) 



HQ) 



Evidently, e*(x,Q)Je(x,Q) = 1. The symmetry 2.4 produces an analog of mon- 
odromy property 2.9 for the function e*: 

e*(x + 2l,Q) = e*(x,Q)JM~ 1 (l,-l,\)J- 1 = w- 1 (Q)e*(x,Q). (2.17) 



7 

'The action of the differential operator D = J2 ujjd° on the row vector /+ is defined as 



f+D= E(-d) j (f + ^). 

3=0 
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Lemma 2.2. The function e*(x, Q) has simple poles at the branch points sj:. It 
has fixed zeros at 71, . . . ,7 9 +i- The other zeros for each component of the vector 
function e* lie on every real oval and depend on the parameter x. The function e* 
has the asymptotics at infinities 

e*(x, Q) = ± e =F9(*+0 [e T Je T + o(l)] , Q G (P + /P_). 

Proof. The function Q(Q) is meromorphic with 2(g + l) poles on both sheets above 
points n n lying in open gaps and 2(g + 1) zeros at the branching points sf. At 
infinities it has the asymptotics 

$(Q) = ±l + o(l), Qe(P ± ). 

Now it is easy to prove properties of the function e*(x,Q). It has poles at the 
branch points which arise from zeros of It has zeros at 71, ... , 7 fl +i, the 

poles of Other poles of $(<5) are annihilated by the poles of e + (x, Q). Other 

g+ 1 zeros of e*(x, Q) which depend on x are produced by the corresponding zeros 
of e + (x, Q). 

The asymptotics follows from the asymptotics for e(x,Q) and <&(Q). □ 
Consider periodic variations of the matrix V(x, A) : V = V + 5V. Then = 

p(Q) + Sp(Q) + . . .. We need a standard formula connecting the variations Sp(Q) 

and SV. 

Lemma 2.3. The following identity holds 

i5p(Q) =< e*(x,Q)J5Ve(x,Q) > . 

Proof. Let e(x, Q) be a Floquet solution corresponding to the deformed potential 
V. From the definition 

e + (x, Q) {[Jd x - JV]e(x, Q)) = 

and 

(e + {x,Q)[Jd x -JV\)e(x,Q) = 0. 
Subtracting one identity from another, we have 

e + (Jd x e) - (e + Jd x )e = e + JVe - e + JVe = e + J5Ve. 

Integrating both sides, we have 

+1 

e + Je' + e + 'je dx = e + Je\t\. 

-1 

Using the identities 

e+(l,Q)=e-^ 2l e + (-l,Q) 

and 

e(l,Q) = e^ 21 e + (-l,Q) 
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for the LHS, we have 

( e ii>(Q)2i e -ip(Q)2i _ ^ Q)je(x, Q) = i5p{Q)2l${Q) + lower order terms. 
The RHS is equal to 

21 < e + (x,Q)J5Ve(x : Q) > +lower order terms. 

Collecting leading terms, we obtain the stated identity. □ 
Consider a real hyperelliptic spectral curve T of finite genus corresponding to 
some periodic potential ip. Let us introduce the Baker- Akhiezer function e(r, x, t, Q) 
which depends on three parameters (times) r, x and t and has the asymptotics at 
infinities 

e(r, x, t, Q) = e ±l H r+ ^ 2 *) x [e /e + o(l)] , Q G (P+/P-). 

The BA function has poles at the points 7's, located on the real ovals. These 
properties define the BA function uniquely. The BA function can be written ex- 
plicitly in terms of theta-functions of the curve T, [10]. The BA function has Bloch 
property in x-variable e(r, /, t, Q) = w(Q)e(r, —I, t, Q) and satisfies the identities 

[Jd T - JVi(t, x, t)} e(r, x, t, Q) = 0, 
\Jd x - JV 2 (r,x,t)]e(r,x,t,Q) = 0, 
[Jd t - JV 3 (t, x, t)\ e(r, x, t, Q) = 0. 

The three matrices Vi, V2 and V3 are given at the beginning of this section. 
Let us define the dual BA function e + (r,x,t,Q) at the point Q as 

e + (r, x, t, Q) = e(r, x, t, e±Q) T . 

The identity w(Q)w(e±Q) = 1 implies e + (r,l,t,Q) = w{Q)- 1 e + (r, -l,t,Q). The 
dual BA function e + (r, x, t, Q) satisfies dual identities e + (r, x, t, Q) [Jd T — JVi] = 
0, etc. 

2.2. Symplectic Structures. We assumed in the previous section that the phase 
space consists of smooth 2/-periodic functions ip(x). Instead we can change the 
language and think about the phase space as a space of operators d x — V2 with 
2/-periodic potential. The flows of the NLS hierarchy act on this space as well 
they act on the space of functions ip. All notions of differential geometry can be 
applied to this space of operators with obvious conventions. 

Lemma 2.4. The formula 

u = ^res < e*J5V A 5e > dX, 

p± 

defines a closed 2-form u on the space of operators d x — V2 with periodic potential. 
The flows e tXm , m = 1, 2, . . . on the space of operators defined by the formula 

[9r m ~ V m , d x - V 2 ] = 0, 
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are Hamiltonian with the symplectic structure u Q and the Hamiltonian function 
H m (up to unessential constant factor) . 

Remark. The formula 

oj n = ^2 res ^ n < e *J$ v A5e> d\, n = 0, 1 . . . , 

p± 

defines a closed 2-form u n on the space of operators d x — V 2 with periodic potential 
which satisfy the constrains H k = const, k = 1, . . . , n; see for details [13]. 

Before proceeding to the proof of the Lemma, we compute the first two sym- 
plectic structures using the formula 

<jj n = ^ res — — — < e + J5V A 5e > dX, n = 1,2. 

The result is 

lu = 2i KStjj ASi/j >, (2.18) 

Ul = Ul =< Sip AS^' + 5^ ASi/j' + 25d~ 1 \ip\ 2 A5\ip\ 2 > . (2.19) 

subject to the constraint Hi = const. We present the computaion divided in small 
steps. 

Step 1. Identity $(r ± Q) = -®{Q) implies 

= <Ao + f + !» + •■■, Qe(P + ), 

Using definition of e(x, Q) and e + (x,Q) from previous section, we have 

00 = 1, 

01 = — < di + di > . 



Step 2. Near P + we obtain 



<e+J5\/A5 e >=^ + § + ^ + ---, 
A X z X 6 

where 

si = < e^JSV A Sei >, 
s 2 = < e^JSV A 5e 2 > + < e[j5V A 8ei > . 
Similar, at P_, we obtain 

<e+7*VA*e >= - f| - £| . 

A X 1 X 6 
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Using the expansion for $(Q) from Step 1, we derive 

^0 = 0o(si + Sl), 

= 0O(S2 + S 2 ) + 0l(Sl + Si). 

iS'tep 5. Computing si, we have 

si =< e^J5V A 5ei >=< 5^ A5b 1 > . 

Using the formula b\ = —iip, we obtain si = —i < 5ip A 5ip > and 2.18. 
Step 4- The first term in the formula for s 2 produces 

< e^JdV A 5e 2 >=< Sip A 5b 2 > . 

Using recurrence relation 2.14: 6 2 = ib[ — iipdi, we obtain 



The constraint H 1 = const implies < 8\ip\ 2 >= and using the explicit formula 
for di, we obtain 2.19. 

Proof. Closeness of the form uio follows from the result of next Lemma or from 
explicit formula 2.18. For the second statement we present a complete proof only 
for m = 0. Higher flows can be treated similarly. In order to prove that the first 
flow is Hamiltonian one has to establish, [1]: 



where ig t is the contraction operator produced by the vector field Xi. Using time- 
dependent BA functions ia t 8e = e', ia t SV = V we have 



Let us compute the residue at P + . From the computation preceding the proof 



< 5ip A i8b[ — i5ipdi — iip5di > . 
The second term in the formula for s 2 produces 

< e\j8V A 8e 1 >=< A 5d 1 + d^ip A 5bi > . 

Finally 

s 2 =< 5ip A 8~4> - i(di + di)5i) A dip - i5\ip\ 2 A 8d\ > . 
Using the formula for s±, 



Ul =< 8i\) A Sij! + 5ip A 5ip' - i5\ip\ 2 A 8{d l - d x ) > . 




res < e* JV'Se > dX - res < e* J8Ve' > dX. 



1 




Q e (P+), 



HQ) 



l 



Q e (P_). 



Then, using V = [f 03, V] = io~ 3 V, we obtain 

res < e* JV'5e > dX = i(f) < e^Ja 3 Y 5ei > . 
Similarly, using e* = \o 3 e we have: 

res < e*J5Ve' > dX = i2<p 1 < elJSY a 3 e > 

+ i20 o [< e^JSYoa^ > + < eJj5Y a 3 e >] . 

The first term vanishes and 

ig t res < e* J5V A 5e > dX = i(f) < e^Ja 3 Y Q 5ei > 

- -0o [< eg J5Y a 3 ei > + < eJj5Y a 3 e >] ■ 

Similarly at P_ 

%Q t res < e* J8V A 5e > d\ = —i(f)o < Ja 3 Y 5ei > 

+ -0o [< elJ5Y Q a 3 ei > + < e\ JSY a 3 e >] . 

Finally, we obtain 

i 9t ^o = «0o [< eg Ja 3 Y 5e 1 > - < e^Ja 3 Y 5ei >] = -5 2H ± . 

□ 

2.3. Darboux coordinates. The formulas 2.18-2.19 give examples of symplec- 
tic forms. All these forms can be put in the Darboux form in the coordinates 
associated with poles of the Baker-Akhiezer function. 

Lemma 2.5. The formula 

f (Q) =< e*J5V A 5e > dX 

defines meromorphic in Q differential 2-form on F with poles at 71, . . . , 7 9 +i and 
P + ,P-. The symplectic 2-form defined by the formula 

uo = J^res £o(<2) 

p± 

can be written as 

2 9+1 

^0 = - ^2 S P(^k) A <JA(7 fc ). 
fc=i 

Remark 1. The meaning of the right-hand side of this formula is the following. 
The curve V (or its cover T) is equipped with two meromorphic functions X(Q) and 
p(Q). Their variation 5p(Q) and 5X(Q) at the points of the divisor is computed 
for variation of the potential ip(x), ip(x); —l<x<l. The RHS of the formula is 
the sum of an exterior products of these variations. 

20 



Remark 2. In fact for a general smooth potential the divisor A( 7 fc) and values of 
the quasimomentum cosh -1 A(A( 7 fc)) with suitably chosen sign (=sheet) determine 
the potential. In other words they are global coordinates on the phase space. First 
note that the discriminant A can be reconstructed from this data using Shannon 
interpolation, see [19]. Thus the curve V is known. The potential can be effectively 
recovered from the divisor via trace formulas, see [18]. 

Proof. Note 8V = SY does not depend on A. Essential singularity of the Floquet 
solutions at P± cancels out, and infinities are simple poles for the form £o(Q)- 111 
the finite part of the curve £o(Q) has two sets of poles. One is the poles 71, ... , 7s+ i 
of the Baker-Akhiezer function. Another is the branch points of the curve V . By 
the Cauchy theorem 

X/ es = -J^res f (Q) - res £ (Q). 



p± 



Sk 



Let us compute contribution of the first set of poles. Near 7^ we have 

res e 



Therefore, 



and 



Se = 



e = 



res e 



A - A( 7fc ) 



+ 0(1). 



(A - A( 7 *)) 



; 5\( lk ) + O 



A - A( 7fc ) 



5e 



-5A( 7fe ) + O 



Note that e*(x, 7^) 



A-A(7 fc ) \\-\(>y k ) 
and from Lemma 2.3 we obtain 



(2.20) 



res 7fc £ (Q) =< e*J5Ve > (j k ) A 5\(j k ) res 



dX 



A - A( 7fc ) 



iSp(-f k ) A 5A(7 fe ) 



Now consider branch points of the curve. They produce a nontrivial contribu- 
tion, though the pole of e* at the branch point is annihilated by the zero of the 
differential dX. Nevertheless, the variation 5e(x,Q) has a simple pole at s k - First, 
let us make a general remark. 

Consider, the variation of a function f(Q,?p,ip) under variation of the potential 
ip(x), ip(x), —l<x<l; taken for Q in the vicinity of the branch point s k and 
a fixed value of A. Such variation will have a pole at the branch point itself. At 
the branch point A fails to be a local parameter, but w is fine due to the fact 
d\R(X, w)\ Sk 7^ 0. Now, consider a function f(Q,ip, ip) = f(w, ip, ip), and define its 
variation 5 for a fixed value of w. Then, 

Sf = $of + j- Sw. 
dw 
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Take, for example f(Q) = X(Q), then 



= 5 \ + 5w. 
dw 



Therefore, for a general / we have 



5f = 6 f- d l6 X = - d l6 X + O(l). 



(2.21) 
(2.22) 



The zero of the differential dX at the branch point produces the pole of 5f. 

We can proceed to the computation of the residues of C,o{Q) at s^. In the local 
parameter (A — A(sfc)) 1 / 2 ~ w — 1, A = X(Q): 

ei(x) SX(s k ) 



5e(x,Q) 



2 (A-A(a fc )) 1/2 



+ . . . 



where 



e(x, Q) = e (x) + ei(x)(A - A(s fc )) 1/2 + . . . , 



Q e (s fc ). 



Similarly, 



and we have 



ei(x) 



dA 



2 (A-A( Sfe ))V2 



+ 



5e(x,Q) = -^^5X(s k )+0(l). 



dX 

The leading term is the same as in general formula 2.22. Therefore, 2.21 implies 
res Sfe £ (<5) = -res Sfc [< e*J5Vde >] A <5A(s fc ) = res Sfe ' 



< e*J5Vde > A 



dw 



Now using 

e*(x,Q) = e*(l,Q)M T (l,x,X)-\ 

de(x,Q) = dM(x,-l,\)e(-l,Q)+M(x,-l,\)de(-l,Q) 
with the help dM(x,y, X)\ Sk = we obtain, 



res Sfc £ (<2) = res 



e*(Z, Q) < M J (/, x, X)- L J5VM(x, -I, A) > de(-l, Q) A 



c?A 5w 
dw 



Symmetry 2.4 of the monodromy matrix implies 



1 



< M T {1, x, X)~ 1 J8VM(x, -I, A) >= —JSM(l, -I, A), 



and using skew-symmetry of the wedge product 
res Sfc £o(Q) 



21 k 



e*(l, Q)J5M(l, -I, X)de(-l, Q) A 



dX 5w 
dw 



1 
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res 



Sk 



e*(l, Q)J(5M(l, -I, A) - 5w)de(-l, Q) A 



c?A 5w 
dw 
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Identities 2.9 and 2.17 imply 

e*J(5M-5w) = 5e*J(w-M), 

J(w - M)de = J(dM - dw)e. 

Therefore, 



res Sfc f (Q) = Yl TeSsk 



Se*(l)J(dM(l, -I, A) - dw)e(-l) A 



dX Sw 
dw 



Since e*(l) Je(— I) = w 1 , we have 



res Sfc ^o(Q) = Yl TeSsk l e *( l ) JSe (- 1 ) A *«"iA] . 

The one form 

e*(l, Q)J5e(—l, Q) A ^(<2) rf A(Q) 

is holomorphic (in the parameter A) outside of the poles 7^ and the branch points. 
At infinity the essential singularity cancels out and due to 2.6-2.7, Lemmas 2.1 and 
2.2 

e*{l,Q)J5e(-l,Q) A 5w(Q) = o Q 

This implies 

res P± [e*(l,Q)JSe(-l,Q) A 6w(Q)d\(Q)\ = 0. 
By the Cauchy theorem, 

J^res Sfc [e*(l)J8e(-l) A d\ 5w] = - ^res 7fc [e*(l)J8e(-l) A d\ 5w] 

Sk 7fc 

Therefore, using e*(7) = w~ l e*(— I), we have 

X/ eSs * £o(Q) = -^]^2 Tes -/k [e*(l)J6e(-l) A dX 5w] 



Ik 



9.7 



2i^ res 7 fe 

7fc 



e*(-Z)J5e(-Z) A dA 



5w 



w 



Using the formula 2.20 we finally obtain 

X/ eSs * = ^^p(7fe) Ares 7fc [e*(-Z)J5e(-Z)dA] 



7fc 



= X] ^(7fc) A <^(7fc) res 

7fc 



dA 



A - A( 7fc ) 



ik 



□ 
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Remark. For the higher symplectic structures an analogous result holds 
2 9+1 

= -J2 An M7fc) A <JA( 7fc ), n = 0, 1 . . . ; 



fc=l 



subject to the constrains H k = const, k — 1, . . . , n. 

2.4. Action-angle variables. Here we describe briefly another system of Dar- 
boux coordinates. We refer to the paper [19] for details. 
The actions Ik, k — 1, . . . , g + 1; are defined by the formula 

1 



p(\)d\. 



In the formula above the multivalued function p(X) is normalized such that p(sk) 
0, k — 1, . . . , g + 1. The angles 9 k , n = 1, . . . , g + 1; are 



n=l J Sn 



The differentials n = 1, . . . ,g + 1; are of the third kind with poles at the 
infinities P± normalized in such a way that 8 



/ 

Ja k 



a n = 2n5t 



As it is proved in [19] by a direct computation 

R, /*} = <& 

All other brackets vanish 

{/„,4} = 0, {9 n ,9 k } = 0. 

These formulas imply the identity for symplectic forms. 

9+1 
n=l 

We conclude this section with a few remarks. 
Remark 1. Another way to prove the identity for symplectic forms without em- 
ploying the Poisson bracket is found by Krichever, [11]. 
Remark 2. McKean, [17], proved various identities for 1-forms. 
Remark 3. In the finite-gap case the curve Y is specified by 2(^ + 1) branch points. 
The actions Ik, k — 1, . . . , g + 1; together with the g other periods (see 2.8) of 
the differential dp and the constant Pq determine the curve. This fact is due to 
Krichever, [12] (see also [2]). In the infinite-gap case it is shown in [19] that J — 0's 



8 <5^ is Kronecker delta. 
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are global coordinates on the phase space for a general square integrable potential. 
This property holds for a finite gap potential as well. 

3. The Scattering Case 

3.1. Jost solutions. In the next sections we consider the scattering problem 
for the Dirac operator on the entire line with rapidly decaying potential. The 
Riemann-Hilbert approach to the scattering theory for canonical systems with 
summable potential was constructed by M.G. Krein and P.E. Melik-Adamian, 
[8, 9, 16]. This approach was used many times in soliton theory, [4]. 

To simplify the estimates we assume that the potential ip is from the Schwartz' 
space S(M) of complex rapidly decreasing infinitely differentiable functions on the 
line such that 

sup|(l + :r 2 ) n ^ (m) (a;)| <°o m, n = 0, 1, . . . . 

X 

Let us introduce the reduced transition matrix T(x,y, A), x > y\ by the formula 
T(x, y, A) = E- 1 (^j M(x, y, A)^ 1 (-*f) , (3.1) 

where -E(4f ) = exp (— ^cr 3 ) is a solution of the free equation (ip = 0). The matrix 
T(x, y, A) solves the equation 

T'(x, y, A) = Y (x)E(\x)T(x, y, A), T(y, y, A) = /. 

The spectral parameter enters multiplicatively into the RHS of the differential 
equation. The solution is given by the formula 

T(x,t/,A) = exp I* Y (£)E()4)dZ. (3.2) 



The symmetry of the matrix Y : a{Y (x)ai = Y (x) is inherited by unimodular 
matrix T: 

a 1 T(x : y,X)a 1 = T(x,y, A). 

For real A the formula 3.2 and the rapid decay of the potential imply an existence 
of the limit 

T(A) = \imT(x, y, A) = ( J , when y — > — oo and x — > +oo; 

and |a(A)| 2 - |6(A)| 2 = 1. When the potential ip e S{M), we have b{\) e S(R). 

We introduce Jost solutions J±(x,\) as a matrix solutions of the differential 
equation 

J'±(x, A) = V(x, A) J±(x, A), J±(x, A) = E ( — j + o(l), when x — > ±oo. 
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An existence and analytic properties of the Jost solutions follow from the integral 
representations 

J + (x, A) = E (?f) + £°° T + (x, $E (f) df, 

J.(x,X)=E^fj + f_J-ME (f) de. 

The kernels r± are unique and infinitely smooth in both variables. Introducing the 
notation J± = j±\ j'^j we see from the integral representations that j_\x, A), 

j+\x, A) are analytic in A in the upper half-plane and continuous up to the bound- 
ary. Also, the columns j-\x,X), j+\x,\) are analytic in the lower half-plane 
and continuous up to the boundary. 

Now we describe analytic properties of the coefficient a(A) of the matrix T(A). 
The monodromy matrix M(x, y, A) can be written in the form 

M(x,y,X) = J+(x)J-\y) = J^x)Jz\y). 

Therefore, 

J- 1 (x)M(x,y,X)J4y) = J+\y)J-{y) = J?(x) J_(x). 

The variables x and y separate and the above expression does not depend on x or 
y at all. By passing to the limit with x — > +00, y — > —00 we have 

T(A) = J-\y) J_(y) = J~\x) J_{x). 

Therefore, 

a(A)=j (1)T (A)Jjf(A). 
The properties of Jost solutions imply that 

• a(A) is analytic in the upper half-plane and continuous up to the boundary; 

• a(A) is root-free; 

• |a(A)| > 1 and |a(A)| 2 - 1 G S(R) for A real, a(A) = 1 + o(l) as |A| — ► 00. 

This coefficient will be used to construct the scattering curve Too. 

Let Poo(A) be such that a(A) = exp(— i2p 00 (X)) for A in the upper half-plane. The 
quantity j9oo(A) in analogous to the quasimomentum studied in the periodic case, 
see [26] 9 . From the properties of a(X) it follow that Poo{X) is analytic in the upper 
half-plane and continuous up to the boundary; %?oo(A) > for OA > 0; Poo(X) = 
o(l) for |A| — > 00; for real A, the density of the measure d/x^A) = Spoo(A) dX 
belongs to S(M). The function Poo(X) can be written in the form 

d^t) 



Pc 



,(A) = i/^ 

71 J t — 



A 



9 The spectral parameter A in the paper [26] has different scaling and should be replaced by j . 
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Expanding the denominator in inverse powers of A, we obtain: 
°° 1 1 r+°° ff jj 



(3.3) 
The 



where H\,Hi and H3 are the integrals introduced above with I = +00. 
expansion has an asymptotic character for A : 5 < argA < 7r — 5, 5 > 0. 

To describe the asymptotic behavior in x of the Jost solutions j+\x,\) and 
j-\x, A) we assume that A is real and fixed. Then, 

(x,X) 



x — > —00 

o(A)/^(x,A)-5(A)A 
/<-(z,A) 



x 



+00 

a(A)/^(x,A) + 6(A)/_ 



(x, A), 



where 



e '2 1 







/_(z,A) 







2 x 



are solutions of the free equation. 

We sketch the derivation of the asymptotics for j + (x, A), when x — > —00. Let 
ip be a potential that vanishes outside the segment [— L, +L\. In this case formula 
3.2 becomes 



(3.4) 



T(A)=exp J Y (Z)E(\£)dS, 
and T(A) is an entire unimodular function of A of the form 

1 ; [ 6(A) a(A) 
From the definition of the matrix M: 

tf\L,\) = M(L,-L,\)rf\-L,\), 

and from 3.1 

" a(A)e- lAL 6(A) 
6(A) a(X)e lXL 

Obviously j { +\L,\) = /_(L,A) and j { +\-L,X) = ci/_ > (— L, A) + c 2 / < _(— L, A) 
with unknown coefficients Ci and c 2 . Formula 3.4 for real A leads to the linear 
system 

Cl 6(A) + c 2 a(A) = 
cia(A) + c 2 6(A) = 1. 

Solving for c's we obtain the stated formula. For a potential with noncompact 
support one has to take L sufficiently large to make the error negligible. 

The Riemann surface is obtained by gluing together along the real line two 
copies of the complex plane (see Figure 2). One copy we call "+" and another 
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M(L,-L,X) = 



Each copy has an infinity P + or P_. The point Q G is determined by A = X(Q) 
and specification of the sheet Q = (A, ±). Let us define for the "+" copy j(x,Q) 
to be j { +\x,X) if 3A > 0; and j {2) (x,X) if 3 A < 0. For the "-" copy we define 

j(x, Q) to be j W {x, A) if 3 A > 0; and j { +\x, A) if 3 A < 0. In the vicinity of P± 
the function j(x,Q) has asymptotics, 



where A = X(Q) and 



Jo 




1 



3o/j + o(l) 



Jo = 



(3.5) 



1 




Therefore, j(x, Q) can be viewed as a BA function for the singular curve T c 
We also introduce the matrix BA function 



H + (X)= \£\X),jf(X) 



and 



H_(X) - 

analytic in the upper/lower half-plane respectively. They are connected by the 
gluing condition 

H-(x, A) = H+(x, X)S(X), where A G R (3.6) 

and the scattering matrix 5(A) 



S(X) = - 



1 b 
-b 1 



The adjoint (dual) Jost solution j + at the point Q is defined by the formula 

j + (x,Q)=j(x,Qf. 

Any Jost solution satisfies [ Jd x — JV] j = 0. By analogy with the periodic case 
one can prove that j + satisfies j + [Jd x — JV] = 0. 



The matrices and H_ are defined as 



J + 



H+(X) = ^H^X) = 



.(2)T 

3- 
J + 



Extending a into the lower half-plane by the formula a* (A) = a(A), we define 

o"3 



HUX) 



a(A) 



HUX) 



for QfA > 0; 



and 



H*_(X) 



CT2, 



a* (A) 



H-(X) 



for 



It is easy to check that the dual gluing condition holds 

H*_(x, A) = S~\X)H* + {x, A), where 



3A < 0. 



A G 



(3.7) 
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and 



s-'W = - 

a* 



1 -b 
b 1 



Next two lemmas state asymptotic properties of Jost solutions which will be 
used in computations with symplectic forms. 



Lemma 3.1. (i) For fixed x the following formulas hold 

oo oo 

j+\x, A) = e^^ is(x) r = e+^Yl 

s=0 s=0 

where go = 0, k = 1, and 10 

oo oo 



9s 



s=0 



s=0 



k s 
9s 



where g = 0, fc — 1- The expansion has an asymptotic character for A 
arg A < 7r — 5, 5 > 0. 

(wj T/ie coefficients g±, k\ are given by the formulas 



5 < 



and 



9i = 



fci 



r+oo 

J X 



\ 2 dx' 



ki = i 



(x')\ 2 dx'. 



Proof (i). Using 11 dgr+\x, £)l£=oo = 0, for n — 0, 1, ... and integrating ra times 
by parts, 



A 2) (*,A) 



e 2 



= e 2 




1 


1 



/•oo 

+ / rf (z,£)e^ 

J X 



I X 

iXx iXx 

e * „(2), | , e a ( 2 ) , 



e 2 



+(-i) n ^r lr i 2) (^^)+(- i ) n -^ 



This implies the existence and asymptotic character of the expansion in the pa- 
rameter A. The other infinity can be treated similarly. 



10 Operation applied to a scalar signifies complex conjugation and reversal of infinities, see 
formulas in the part (ii) below. 

n r = [r« r( 2 )]. 
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(it). Consider j+\x, A) first. The differential equation j' = Vj implies 



9s 
k' 



+ Y 



9s 
k s 



(I + 0-3) 



9s+i 
k s +i 



s = 0,l,...; 



and g = 0, k = 1. 

This recurrent relation leads to the identities 



Ws+l, 

0. 



-g's + #« 

~K + ipg s 

For s > 1 we have the boundary condition 

9s (x) I x=+oo (^) I x=+oo 

These imply the stated formulas for g±, k\. 
For j ( L > (x, A) the differential equation implies 



k' 



9 



+ Y 



k s 
9s 



2 (- J + *s) 



#5+1 



0,1, 



and fc = 1, go — 0. 

The recurrent relation produces the identities 

-g's + = -ig 8 +i, 
-K + fas = 0. 
For s > 1 we have the boundary condition 



& 0*0 Is 



k s (x) | a 



0. 



These imply the stated formulas for gi, fci. We are done. 



□ 



Remark. It is interesting to compare asymptotic expansions for and 
e(x, Q). For the Jost solution j^l\x, A) normalized at the left 



iXx 
2 



1 





+ 



A 



^■0 



+ 



and 



1 





1 

+ A 



-iipo + i f_ x |^| 2 



+ 



Qe(P- 



If ^ is compactly supported and I becomes sufficiently large, then e t ^ l e(x, Q) = 
j(x,X). For the Jost solution normalized at the right the situation is slightly 
different. If one defines the new e(x,Q) = e(x, Q)w~ 1 (Q) which is the Floquet 
solution normalized at the right end x — I of the interval, then for compactly 



supported ip and sufficiently large / the new e~ l ^ l e(x, Q) = j+\x, A). 

A result similar to Lemma 4.1 holds for Jost solutions analytic in the lower 
half-plane. 
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Lemma 3.2. (i) For fixed x the following formulas hold 

oo oo 

j+\x, A) = e-^j^j.^A- = e-^J] 

s=0 s=0 

where h = 1, / = 0, and 

oo oo 

A) = e**^ j,(s)A- = e + ^^ 



A" 



s=0 



fs 



A" 



where h — 1, / = 0. TTie expansion has an asymptotic character for A 
-5 > argA > -7r + 5, 5 > 0. 

(zzj TTie coefficients hi, f\ are given by the formulas 

"+00 

A = #, h = -i I Mx')\ 2 dx' 



and 



fi = 



"'I 



Similar to the periodic case we will need time dependent BA functions. They 
are obtained by an elementary construction. 

Lemma 3.3. [7]. There exists the Jost solution j(r,x,t,Q) on the curve 
with three time parameters r, x and t which satisfies the differential equations: 

[d r -V 1 (r,x,t)]j(r,x,t,Q) = 0, 

[d x ~ V 2 (t, x, t)]j(r, x, t, Q) = 0, 
[d t -V 3 {r,x,t)]j{r,x,t,Q)=0. 

.(2) 

Proof. We consider "+" sheet and the upper half-plane where j(x, Q) = j + (x, A). 
First, we construct j(t, x, t, Q) such that [d x — V^r, x, t)] j(r, x, t, Q) — normal- 
ized for all r and t as 



j(T,x,t,Q)~e +l ^ 




1 



when x — > oo. 



Then, we construct j(r,x,t,e±Q) on the lower sheet as a solution 

[d x - V 2 (t, x, t)] j(r, x, t, e±Q) = 
for the same value of the spectral parameter A = X(Q) normalized for all r and t 



as 



j(T,x,t,e±Q) ~ e ! 2 : 



1 
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when x — > — oo. 



The solutions j(r, x, t, Q) and j(r, x, t, e±Q) span the kernel of the operator [d x — 
V 2 (r, x,t)]. Now we introduce 

1 A 2 

3new( r , X, t, Q) = e _ *3 T -*T*j( r> Xj t, Q) , 

which is the desired solution. Evidently, 

[d x - V 2 (r, x, t)} j new (r, x, t, Q) = 0. 

To prove the first identity of the statement we note, that commutativity of the 
operators d T — V 1 and d x — V 2 implies 

[d T ~ V±(T,x,t)\ j new (r,x,t,Q) = c 1 (r,t)j(r,x,t,Q) + c 2 (r,t)j(r,x,t,e ± Q). 

From another side as x — > +00, 



[d T - Vi(t, x, t)] j new (r, x, t, Q) = 



a 1 



e -i±T+i^X-i^-t 




1 



+ =0(1). 



Due to the linear independence of the solutions j(r,x,t,Q) and j(r,x,t,e±Q) we 
have ci(r,t) = c 2 (r,t) = 0. 

Similarly it can be proved that 

[d t - V 3 (t, x, t)} j new (r, x, t, Q) = 0. 

Another sheet of can be treated the same way. We are done. □ 
Remark. It is easy to see that, 

j new (T, X, t, Q) = e ±*(-^+^^) \j o /j Q + (1)] Q E ( P± ). 

Thus the standard Jost solution with asymtotics 3.5 can be obtained from the BA 
functions if one puts r and t equal to 0. 

3.2. The symplectic structures. We are ready to introduce the scattering ver- 
sion of the Krichever-Phong formula. The everaging is defined now as an integral 

over the entire line < • >= f + °° dx. 

j —00 

Theorem 3.4. The formula 

u = trace res ]- [< H* + J5V A 5H + > + < H1J5V A 5H_ >] dX, (3.8) 

defines a closed 2-form uo$ on the space of operators d x — V 2 with potential from 
the Schwartz class S(M). The flows e tXm , m = 1, 2, . . . on the space of operators 
defined by the formula 

[d Tm -V m ,d x -V 2 ]=0, 

are Hamiltonian with respect to the 2-form u with Hamiltonian function H m (up 
to a non-essential constant factor). 
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Remark 1. The symbol 

res < H* + J5V A 5H + > 

means the coefficient corresponding to the term j in the power series expansion 
near infinity in the upper half-plane. The second term 

res < H*_J5V A 5H_ > 

is defined in the same way, only the upper half-plane plane is replaced with the 
lower half-plane. 
Remark 2. The formula 

\n 

uo n = trace res — [< H* + J5V A 5H + > + < H*_J5V A 5H_ >] dX, 

where n — 0,1 ... , defines a closed 2-forms u n on the space of operators d x — V2 
with potential from the Schwartz class S(OL) that satisfy the constraints Hk = 
const, k — 1, . . . ,n. It is instructive to compute explicitly the symplectic forms 
for small n. The first few are given by the formulas 

ujq = 2i < 5ip A 5ip >, 

and 



ui =< 5ip A 5ip' + 5ip A 5ip' + 5 



01-/ 101 



A^l^l 2 > 



subject to the constrain Hi = const. The derivation employs Lemmas 3.1-3.2 and 
similar to the periodic case. 

Proof. Closeness of the form u>o follows either from the explicit formula or from 
the result of the next theorem. The proof of the second statement we present 
for the first e tXl flow. The time dependent Jost solutions entering into 3.9 are 
constructed in Lemma 3.4. Let ig t be the construction operator produced by the 
vector field Xi. We will prove ig t ouo = —52Hi. For the first term in 3.8 

trace res < H* + J5V A 5H + > dX 

= res i < j ( ^ )T J5V A Sj { + ] >dX-res^< jf T J5V A 5j W > dX[3.9) 
Applying the contraction operator to the first term in 3.9 



i dt res - < j m J5V A > dX 

CL 



res - < j WT JV5£ ] >dX-res -< £ )T JSv£ ] ' > dX. 

a a 



From 3.3 we have 



1 a,i _ 

- — a + — + ... , where a = 1, a x = — 1 

a X 



too 

L/.I2 



—00 
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Using V* = [|cr 3 , V] = ia 3 V ', and Lemma 3.1, we have 

res - < j (1)T 'jV5j { + >d\ = a < jljicr 3 V5j l >= - < if}8if> > . 

(X 

Similarly, using j' = ^cr 3 j, we have 
res I < jWjSVj®' >d\ = ai < j T JSV l -a 3 j > 



+ a 



a T % ^ T l 

< Jo J5V-a 3 j 1 > + < j 1 JSV-a 3 j > 



The first term vanishes, the second produces 
Finally, 



- < if}5if} — ipSif} > . 



i dt res - < j a)T J5V A 5jf >d\ = -5Hl 

(X 

The second term in formula 3.9 can be treated similarly. Therefore, 

trace res < H* + J5V A 5H + >= -25Hi. 
The second term in formula 3.8 produces the same result. The proof is finished. □ 

3.3. Action-angle variables. In response to infinitesimal deformations of the 
matrix V = V + SV the matrix T(A) changes according to the rule: T(A) = 
T(A) + ST(X) + . . .. The next result is similar to Lemma 2.3 of the periodic case. 

Lemma 3.5. The following formula holds 

—8a b5a — ab% 



< H+J5VH+ > 
with averaging defined as 



a5b — b5a —5a 



/+oo 
H+(x, X)J5V(x)H + (x, A) dx. 
-oo 

Proof. Let us assume, first, that if> has compact support. We denote by V, T and 
j deformed matrices V,T and the Jost solution j. We will derive the expression 

• (2) T '(1) 

for < j\' J5Vj + > in the left-upper corner. First, we obtain the formula 

_ +L 

< j + J5Vj > + lower order terms = j + Jj . (3.10) 

— L 

Indeed, 

j + {[Jd x - JV]j) = 0, 
(j + [Jd x -JV])j = 0. 
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Subtracting one identity from another we have 

j + JSVj + lower order terms = j + (jdj^j — (j + Jd) j. 
Integrating the RHS in x variable we obtain 



+ L r n _ +L 

3 + Jj + 3 + 'Jj\ dx = j + Jj ^ 



This implies 3.10. Now using formulas for the asymptotics of j + and j for the 
RHS of 3.10, we have 



~af^ + bf. 



- [afl - bfl] J/. 



af^Jf^ — afi^.Jf^_ = a — a = —5a + lower order terms. 

Collecting terms of the same order, we obtain the result. The case of a potential 
with non-compact support can be considered using approximation arguments. For 
other entries the arguments are the same. Lemma is proved. □ 

Theorem 3.6. The following formulas hold 

+oo _ 



J_ r 5b(X) 
iri J \a( 



ni J | a( A) j 

— oo 

Proof. By the Cauchy integral formula 

+oo 

1 1 f 

- trace res < H1J5V A 6H+ >= -— / trace < H1J5V A 5H + > dX, 

2 2m J 

— oo 

and 

+oo 

- trace res < H* J5V A5H_>=— / i racr < // JA ' A W/ > <1\. 

2 27ri 



- / 

>7T2 J 

Taking sum 



1 



uo = / trace < H*_J5V A 5H_ > d\ / trace < H*J5V A 5H+ > dX. 

2ixi J 2iri 1 

— oo 

Using 3.6, 3.7 

8H- = 5H+S + H + 5S, 
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we obtain 



-t-oo i-oo 

u = ^- J trace S' 1 < H* + J5V A5H + > S d\+ J trace S' 1 < H* + J5VH + > A5S dX 



-oo 

+00 



1 



2m 

—00 

+00 



trace < H* + J5V A 5H+ > dX 



J trace < H* + J5VH + > ASSS^dX. 



Now, applying the result of Lemma 3.5, we have 

A 8S S -1 d\. 



+00 

uj = - — ; / trace 



2m 



—5a b5a — a5b 
a5b — b5a —5a 



—00 



After simple algebra we arrive at the stated identity. Theorem is proved. □ 
Remark 1. The formula of Theorem can be put easily into more familiar form 
using the identities 



a| 2 -|6| 2 = l, 51og|6| 



2| 


a\8\a\ 


a 


I 2 - 1 



Indeed, 



1 5b(X) A 5b(X) |&(A)| 2 ^log |&(A)| 2 A ^phb(A) 
7 |a(A)| 2 ~ |a(A)| 2 

= 251og|a(A)| A<5phfe(A). 

Therefore, 

+ OO _ +00 

1 f 5b(X)A5b(X) 2 f fl , , x . ( . , L/x . M 

^o = — / 1 -jT^r = ~ / ^ log a( A) A 5 ph 6( A) dX. 
m J \a(X)\ 2 tc J 

—00 —00 

Remark 2. The formula 

+00 _ 
1 f 5b(X) A 5b(X) 
u n = — I / ■> \ 12 ^ dX, n=l,2,...; 

m J \a(X)\ 2 

—00 

subject to the constrains Hk = const, k = 1, . . . , n; gives Darboux coordinates for 
higher symplectic forms. 
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